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INTRODUCTION AND IMPORTANCE OF HIGH-REDSHIFT X-RAY STUDIES 



ABSTRACT 

>^ I X-ray studies of high-redshift {z > 4) active galaxies have advanced substantially over the past few years, 
largely due to results from the new generation of X-ray observatories. As of this writing X-ray emission has 
• been detected from nearly 60 high-redshift active galaxies. This paper reviews the observational results and 
' their implications for models of the first massive black holes, and it discusses future prospects for the field. 
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2 Understanding of the X-ray emission from z > 4 active galactic nuclei (AGN) has advanced rapidly over 
I the past few years. The high sensitivities of Chandra and XMM-Newton have allowed the efficient detection 
' of many z > 4 AGN, and wide-field AGN surveys (e.g., the Sloan Digital Sky Survey, hereafter SDSS; York 
!> ■ et al. 2000) have greatly enlarged the number of suitable X-ray targets. The number of X-ray detections 
^ . at z > 4 has increased to 57 (see Figure la).' This increase has allowed the first reliable X-ray population 
^ I studies of z > 4 AGN, and it has been possible to obtain X-ray detections at redshifts up to z = 6.28 (e.g., 
Brandt et al. 2002a). 

The X-ray emission from AGN provides direct information about their black hole regions, where accretion 
and black hole growth occur; the emission is thought to be produced by the inner accretion disk and its 
corona (e.g., Poutanen 1998). Do the early black holes seen at high redshift feed and grow in the same way as 
local ones? X-ray observations of high-redshift AGN can address this fundamental question. It is plausible 
that high-redshift AGN could be feeding and growing differently. The comoving number density of luminous 
quasars changes by a factor of > 100 over the history of the Universe (e.g., Figure 9 of Fan et al. 2001 
and references therein), and part of this strong evolution is believed to be due to environmental changes 
that might also impact the X-ray emission region. There have also been theoretical predictions that AGN 
accretion rates, relative to the Eddington rate, should change with redshift (e.g., Kauffmann and Haehnelt 
2000); such changes can lead to accretion-disk instabilities (e.g., Gammie 1998 and references therein) and 
radiation "trapping" effects (e.g., Begelman 1979). 

X-ray absorption measurements can also be used to probe the large-scale environments of high-redshift 
AGN. Changes in the amount of X-ray absorption with redshift have been discussed by many authors. 
The fraction of radio-loud quasars (RLQs) with heavy X-ray absorption appears to rise with redshift, with 
column densities of >2 x 10^^ cm^^ being seen at 2;>3 (e.g., Elvis et al. 1998; Fiore et al. 1998; Reeves 
and Turner 2000). The absorbing gas may be circumnuclear, located in the host galaxy, or entrained by the 
radio jets. Radio-quiet quasars (RQQs) definitely show less of an absorption increase with redshift than do 

'See ,http://www.astro.psu.edu/users/niel/papers/highz-xray-detected.dat, for a regularly updated list of z > 4 X-ray detections. 
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Fig. 1. (a) ^i?i45o(i+2) magnitude versus redshift for z > 4 AGN (small solid circles). Solid triangles (large solid 
circles) indicate AGN that have been detected by Chandra or XMM-Newton {ROSAT). Open triangles (large 
open circles) indicate AGN that have tight X-ray upper limits from Chandra or XMM-Newton [ROSAT). (b) 
Observed-frame, Galactic absorption-corrected 0.5-2 keV flux versus ^-61450(14-2) magnitude for 2: > 4 AGNs. 
Object types are shown in the key; see Figure 4 of Vignali et al. (2003a) for further explanation. The slanted 
lines show aox=— 1-5 and aox=— 1-8 loci at a fiducial redshift of z = 4.6. 



RLQs (e.g., Fiore et al. 1998). However, as discussed below, the X-ray absorption properties of liigh-redshift 
RQQs are only now becoming clear. 

We adopt Hq = 65 km s~^ Mpc~^, Qm = 1/3, and = 2/3 throughout. This review partially updates 
Brandt et al. (2002b); some of the technical details and references in that review are not repeated here. 

RECENT X-RAY OBSERVATIONS AND RESULTS 
X-ray Observations 

Most z > 4 X-ray detections to date have been obtained via moderate-depth [~ (1-20) x 10~^^ erg cm~^ s~^ 
in the 0.5-2 keV band] follow-up observations of AGN discovered at other wavelengths (e.g., Kaspi et al. 
2000; Vignah et al. 2001; Brandt et al. 2002a; Bechtold et al. 2003; Vignah et al. 2003a, 2003c). ^ Due 
to its superb imaging capability, Chandra is particularly effective at such observations; a measurement of 
only ~ 2-3 photons at the known position of a quasar constitutes a significant Chandra detection! Chandra 
"snapshot" observations of only 4-10 ks are now reliably yielding X-ray detections of many z ~ 4-6 quasars. 
For example, we have made 30 such snapshots of 2: > 4 quasars and have detected 87% of our targets; we 
have 25 more snapshots scheduled for Chandra Cycle 4. Our targets have included the highest redshift 
SDSS quasars (2^4.8; Mb ~ —27.0 to —28.5) and the most luminous Palomar Digital Sky Survey (DPOSS; 
Djorgovski et al. 1998) quasars (z ~ 4.1-4.5; Mb ~ —28.5 to —30.5). We are also now targeting RLQs 
and "exotic" quasars (e.g., the quasars lacking emission lines in Fan et al. 1999 and Anderson et al. 2001) 
at 2; > 4. Whenever possible, we have obtained near-simultaneous optical imaging and spectroscopy with 
the 8-m Hobby-Eberly Telescope (Ramsey et al. 1998), so that the X-ray and optical properties can be 
compared without concerns about variability. 

Seven 2 > 4 AGN have been selected via their X-ray emission, spanning the redshift range z = 4.14-5.19 

^In a handful of cases, moderate-depth X-ray data were serendipitously available for a field in which a z > 4 quasar was 
subsequently discovered (predominantly using data at other wavelengths). 
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(Henry et al. 1994; Zickgraf et al. 1997; Schneider et al. 1998; Barger et al. 2002; Silverman et al. 2002; 
Castander et al. 2003). These objects suffer from less selection bias than AGN found at other wavelengths; 
the 2-40 keV X-rays detected from z > 4 AGN are extremely penetrating and thus unlikely to be absorbed. 
In addition, the z > 4 X-ray selected AGN found in deep X-ray surveys (i.e., the Chandra Deep Field-North 
and the Lockman Hole) arc substantially less luminous (with Mb ~ —21.4 to —24.4) than those found by 
the SDSS and DPOSS. These moderate-luminosity AGN are likely to be the majority AGN population at 
z > 4; as such, they deserve at least as much attention as the highly luminous SDSS and DPOSS quasars. 

Basic X-ray Properties: Fluxes and Luminosities 

Figure lb shows the X-ray and optical fluxes for z > 4 AGN with sensitive X-ray observations. Most 
of these objects, including even the optically brightest ones known, arc faint X-ray sources with 0.5-2 keV 
fluxes of < 3 X 10^"'^'^ erg cm^^ s^^. The only exception to this itiIc is a small group of radio-loud blazars 
which have 0.5-2 keV fluxes up to ~ 6 x 10^^'^ erg cm^^ s^^ (e.g., Moran and Helfand 1997; Yuan et al. 2000; 
Fabian et al. 2001a,b). It is also clear from Figure lb that, as expected from low-redshift AGN, there is a 
correlation between X-ray and optical flux; Vignali et al. (2003a, 2003c) find this correlation to be significant 
at the > 99% confidence level and give the best-fit relationship (see their §3.1). The scatter around the 
best-fit relationship, however, is substantial (a factor of w 3 in X-ray flux). Practically, this scatter makes 
it risky to attempt long X-ray spectroscopic observations before an X-ray flux has been established via a 
shorter observation. 

The 2-10 keV rest-frame luminosities of the AGN shown in Figure lb span the wide range ~ 2 x 10^^ erg s~^ 
(e.g., Vignali et al. 2002) to « 2 x 10"^^ erg s~^ (e.g., Fabian et al. 1998). Only the radio- loud blazars inhabit 
the highest X-ray luminosity range; RQQs generally have 2-10 keV X-ray luminosities of < 3 x 10^^ erg s~^ 
(e.g., Vignali et al. 2003a). These luminosities assume isotropic emission and no amplification due to 
gravitational lensing (e.g., Wyithe and Loeb 2002). 

X-ray Contribution to the Spectral Energy Distribution 

In AGN studies the contribution of X-rays to the spectral energy distribution is often parameterized with 
aox, the slope of a nominal power law between rest-frame 2500 A and 2 keV [oox = 0.384 log(/2 kev//25oo a) 
where /a keV is the flux density at 2 keV and /2500 A density at 2500 A].^ AGN with larger 

negative ckox values have relatively weaker X-ray emission. The optically selected RQQs at z > 4 have more 
negative aox values ((oox) = —1.74 it 0.02) than low-redshift comparison samples (e.g., the Bright Quasar 
Survey RQQs have (aox) = —1.56 ± 0.02; Brandt et al. 2000; see Figure 2a). However, the z > 4 RQQs are 
also more luminous on average (at 2500 A; see Figure 2b), and thus one must address if the difference in 
CKox is due to redshift or luminosity (or both). Previous studies of this topic have generally concluded that 
aox primarily depends upon luminosity rather than redshift (e.g., Avni et al. 1995 and references therein), 
although Bechtold et al. (2003) recently argued that aox depends primarily upon redshift. Most work on this 
topic has been limited by heterogeneous AGN selection criteria or limited coverage of redshift or luminosity. 

We have been investigating the aox versus redshift /luminosity issue with a sample of 137 SDSS RQQs at 
z = 0.16-6.28 that lie in pointed ROSAT, Chandra, and XMM-Newton observations (Vignali et al. 2003b, 
hereafter VBS; also see Vignali et al. 2003c). We have focused on RQQs since they comprise the bulk 
(r^ 85-90%; e.g.. Stern et al. 2000) of the quasar population and also are "simpler" systems to study (since 
it is not necessary to decompose the X-ray emission into accretion-disk-linked and jet- linked components). 
Compared to much of the previous work, our analysis has several significant advantages: 

1. The SDSS quasars utilized have been optically selected in a well-defined manner and span large ranges 
of redshift and luminosity (e.g., Schneider et al. 2002 and references therein). 

2. Sensitive radio data from the FIRST and NVSS surveys allow effective selection of only RQQs. 

3. We have removed z > 1.5 Broad Absorption Line quasars (BALQs) and quantitatively constrained the 
effects that unidentified z<1.5 BALQs can have upon our results. The underlying X-ray emission 

^For z > A AGN we have generally calculated /2 keV from the observed-frame 0.5-2 keV flux (assuming an X-ray power-law 
photon index of F = 2). Thus, the derived aox values are actually based on the relative amount of X-ray flux in the 0.5(1-1-0) keV 
to 2(1 + 2) keV rest-frame band. 
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Fig. 2. The parameter ctox versus (a) redshift and (b) 2500 A luminosity for optically selected RQQs. The 
symbols at z > 4 are the same as for Figure lb. The open squares are for seven luminous, absorbed Bright 
Quasar Survey (BQS) RQQs (BALQs and mini-BALQs), and the solid squares are for the other 46 luminous BQS 
RQQs (from Brandt et al. 2000). The solid circles with error bars show stacking results for Large Bright Quasar 
Survey (LBQS) RQQs from Figure 6d of Green et al. (1995). The data points in (a) from z = 2-4 are taken 
from the ROSAT study of SDSS RQQs by Vignali et al. (2003b). 

properties of BALQs are often difficult to measure because their X-ray spectra are "corrupted" by 
heavy absorption (e.g., Green et al. 2001; Gallagher et al. 2002; see the open squares in Figure 2). 

4. We have examined the possible effects of flux amplification due to gravitational lensing of z > 4 quasars 
(e.g., Wyithe and Loeb 2002). aox will probably not be changed by gravitational lensing, but the true 
luminosity will differ from that assumed. 

5. We have adopted the observationally favored "ACDM" cosmology and have assessed the sensitivity of 
our main results to the adopted cosmology. 

Partial correlation analyses, using the method of Akritas and Siebert (1996), find an anti-correlation between 
Oox and 2500 A luminosity at the 3.4-3.9o" level (controlling for the effect of redshift). By contrast, no highly 
significant redshift dependence of Oox is found (1.7-2.3(t; controlling for the effect of luminosity) . Thus, VBS 
find results that are qualitatively consistent with those of most previous workers: Oox appears to depend 
primarily upon luminosity rather than redshift.^ 

For z > 4 AGN, the lack of any strong changes in Oqx with redshift (after controlling for luminosity 
effects) is generally consistent with the lack of strong spectral evolution at other wavelengths. For example, 
the ultraviolet rest-frame spectra of z > 4 quasars are only subtly different from those at lower redshifts 
(e.g., Schneider et al. 1989; Constantin et al. 2002), and the fraction of RLQs does not appear to change 
significantly with redshift (e.g.. Stern et al. 2000). 

X-ray Spectroscopy 

X-ray spectral analyses of z > 4 AGN can determine the photon index of the X-ray power law and the 
properties (column density and ionization level) of any intrinsic absorption. Aside from the small group of 

*We note, however, that there is substantial scatter in Ofox (~ ±0.2) at all luminosities and redshifts; the correlations are far 
from being "tight." 
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Fig. 3. (a) Stacked X-ray spectrum constructed from nine luminous DPOSS quasars at z = 4.09-4.51 fitted 
with a power-law model and Galactic absorption (see the text). Data-to-model residuals, in units of cr, are shown 
in the bottom panel. The insert shows the 68, 90, and 99% confidence regions for photon index and intrinsic 
column density obtained via fitting with the Cash (1979) statistic, (b) Hard X-ray photon index versus redshift. 
The solid circles are the 2-10 keV data points with the smallest statistical errors from Page et al. (2002). 
The stars are, from left to right, the data points from ARM 08279+5255 at z = 3.91 (Chartas et al. 2002), 
nine luminous DPOSS quasars at z = 4.09-4.51 (Vignali et al. 2003a; plotted at their average redshift), and 
CXOHDFN J123647.9+620941 at z = 5.186 (Vignali et al. 2002). 

radio-loud blazars mentioned above, detailed X-ray spectral analyses of individual z > 4 AGN have not yet 
been performed due to their low X-ray fluxes (although see Chartas et al. 2002 and Hasinger et al. 2002 for 
X-ray spectral analyses of the z = 3.91 gravitationally lensed BALQ APM 08279-1-5255). 

It has been possible, however, to derive average X-ray spectral constraints via the joint fitting of several 
z > 4 AGN (using the Cash 1979 statistic). Such joint fitting can provide a more stable estimate of average 
X-ray spectral properties than a single-object spectrum of the same statistical quality (since any single 
object might be atypical). Vignali et al. (2003a) have recently performed joint Chandra spectral fitting of 
nine luminous DPOSS quasars at z = 4.09-4.51; together these quasars have ~ 340 counts, and only ~ 2 
background counts are expected. The results are shown in Figure 3a and constrain the rest-frame 2-30 keV 
emission. The spectra are well fit by a power-law model with only Galactic absorption. The derived power- 
law photon index is T = 1.98 ± 0.16 (errors are for 90% confidence), entirely consistent with those found 
for z ?a 0-3 samples of RQQs in the 2-10 keV band (see Figure 3b; e.g., George et al. 2000; Reeves and 
Turner 2000; Page et al. 2002).^ The average constraint on neutral intrinsic absorption, assuming solar 
abundances, is A^h < 8.8 x 10^^ cm~^. If X-ray absorption were present at the level found in high-redshift 
RLQs, it should have been detected. Vignali et al. (2003c) have also used joint Chandra spectral fitting to 
constrain the average spectral properties of 13 SDSS quasars at z = 4.81-6.28. These objects are well fit by 
a r = 1.84 lb 0.31 power law and Galactic absorption 

We note that Bechtold et al. (2003) have reported X-ray spectral constraints for 16 luminous quasars 
at z = 3.67-6.28. They do not report the results of a joint spectral analysis, but the mean of the photon 
indices in their Table 2 is (P) = 1.50 ± 0.15. This mean photon index is significantly harder (at the ~ 2.5fT 
level) than those reported above for the DPOSS and SDSS quasars. 

Vignali et al. (2002) have performed X-ray spectral analyses of the three z > 4 moderate-luminosity 
AGN found in the Chandra Deep Field-North; such analyses are only possible due to the extremely deep 
2 Ms exposure. The highest quality X-ray spectrum was obtained for the low-luminosity quasar CXO- 

^We note that, at present, there is no evidence that the photon index of the underlying X-ray power law has a luminosity 
dependence (after allowing for Compton- reflection effects). 
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HDFN J123647.9+620941 {z = 5.186; Mb = -23.4; 138 counts) and is well fit by a T = 1.81 ± 0.31 
power-law model with Galactic absorption. The other two AGN, CXOHDFN J123642.0+621331 {z = 4.424; 
Mb = -21.6; 45 counts) and CXOHDFN J123719.0+621026 {z = 4.137; Mb = -21.4; 117 counts), have 
Seyfert-like luminosities. They have flatter X-ray spectra with effective photon indices of F ~ 1.6 and 
F ~ 1.1, respectively, when only Galactic absorption is assumed in the fitting (these fits arc statistically 
acceptable). The X-ray spectrum of CXOHDFN J123719. 0+621026 can also be acceptably modeled with a 
F 2.0 power law and an intrinsic column density of 2 x 10^^ cm~^. 

X-ray Variability 

At present, information on the X-ray variability of z > 4 AGN is very limited. The snapshot X-ray 

observations used to detect most z > 4 quasars to date (sec above) cover timescales of only ~ 20 minutes 
in the rest frame, and most luminous quasars at low rcdshift do not vary strongly on such short timescales. 
Highly significant X-ray variability at z > 4 has only been detected from the blazars PMN J0525-3343 
(Fabian et al. 2001b), GB 1428+4217 (Fabian et al. 1999; Boiler et al. 2000), and GB 1508+5714 (Moran 
and Helfand 1997). CXOHDFN J123719.0+621026 in the Chandra Deep Field-North may also be X-ray 
variable (Vignali et al. 2002). There is a tentative claim that quasar X-ray variability increases with redshift 
in the sense that quasars of the same X-ray luminosity are more variable at z > 2 (Manners et al. 2002). 

X-ray Imaging 

X-ray imaging of ^ > 4 AGN can allow searches for spatial extent due to gravitational lensing, small-scale 
X-ray jets, or X-ray scattering by dust in the intergalactic medium. The only z > A AGN with published 
evidence for X-ray spatial extent is SDSS J103027. 10+052455.0 at z = 6.28, and Schwartz (2002) suggests 
that this extent may be due to gravitational lensing. The putative X-ray extent needs to be confirmed 
with a deeper Chandra observation, and lensing will be tightly constrained with an HST observation. Telis 
et al. (2003) have set a direct upper limit on the density of cosmological dust from the absence of an X-ray 
scattering halo around the z = 4.30 blazar GB 1508+5714; they find that intergalactic extinction is unlikely 
to have significantly affected the Hubble diagram of Type la supernovae. 

X-ray imaging can also search for companion objects such as clustered AGN in large-scale structures (e.g., 
Djorgovski 1999; Djorgovski et al. 2003) or large-scale X-ray jets. Schwartz (2002) tentatively reported the 
detection of a large-scale X-ray jet associated with the z = 5.99 quasar SDSS J130608. 26+035626. 3, although 
such a detection would have been somewhat surprising given that this object (and its putative jet) is not 
a powerful radio source (i.e., its Kellermann et al. 1989 radio-loudness parameter is R < 8.7). Deep 
7-band imaging of the putative jet (Ivanov 2002) suggests the X-ray emission is probably associated with an 
unrelated foreground galaxy. Nevertheless, jets should be searched for from other z > 4 AGN. Vignali et al. 
(2003a) and Vignali et al. (2003c) have searched for statistical overdensities of companion X-ray sources 
near DPOSS and SDSS quasars, but no overdensities are found. 

CONCLUSIONS AND FUTURE PROSPECTS 

From a general perspective, the main conclusions from the X-ray analyses reviewed above are the following: 

1. To the limits of current observation, AGN at 2; « 4-6 and z « 0-3 have reasonably similar X-ray and 
broad-band spectra (after controlling for luminosity effects upon the broad-band spectra). Specifically, 
the distributions of X-ray power-law photon index and aox appear consistent. 

2. The current X-ray data do not provide any hints for different accretion mechanisms at low and high 
redshift. The X-ray data are consistent with the idea that early black holes at high redshift feed and 

grow in the same way as local ones. If the accretion mechanism docs change at high rcdshift, it does 
so without dramatically affecting the X-ray emission made in the immediate vicinity of the black hole. 

3. The small-scale X-ray emission regions of AGN appear relatively insensitive to the strong large-scale 
environmental changes that have occurred from z ~ 0-6. The quasar population evolves strongly over 
this time, but individual quasar X-ray emission regions appear to evolve much less. 
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4. The X-ray data provide support that, at least to first order, the spectral energy distributions of local 
AGN can be used when (1) estimating the black hole masses of high-redshift AGN via the Eddington 
argument [~ (1-5) x 10^ Mq for luminous SDSS quasars; e.g.. Fan et al. 2001] and (2) calculating the 
effects of early AGN upon the intergalactic medium (e.g., Venkatesan et al. 2001; Gen 2003; Machacek 
et al. 2003). 

5. X-ray emission appears to be a universal property of luminous AGN, even at z ~ 4-6. The X-rays 
also do not appear to be generally blocked from view by material in nascent AGN host galaxies.^ 

There are a number of ways the work described above can be extended including the following: 

1. At present, there is only a limited number of X-ray detections at z > 5 (see Figure la). It is important 
to increase this number as further z > 5 AGN are discovered, so that strong statistical conclusions 
can be drawn about X-ray emission at the highest redshifts. For example, the SDSS should find ~ 100 
quasars with z > 5 and ~ 30 with z > 6; it can find quasars up to z ~ 6.5 (e.g.. Fan et al. 2003). 
Hopefully X-ray (e.g., Alexander et al. 2001; Barger et al. 2003) and infrared (e.g.. Warren 2002) 
selection can be exploited to discover even higher redshift AGN. Ultimately, with missions such as 
Generation-X , it should be possible to detect ~ 10^ Mq black holes in proto-AGN at z ~ 10; these 
objects will have X-ray fluxes of ~ 5 x 10^^^ erg cm^^ s~^ (compare with Figure Ib).'^ 

2. Minority AGN populations at z > 4 need to be investigated better in the X-ray regime. These include 
quasars lacking emission lines (e.g.. Fan et al. 1999; Anderson et al. 2001; Vignali et al. 2001), BALQs 
(e.g., Brandt et al. 2001; Goodrich et al. 2001; Vignah et al. 2001, 2003a, 2003c), and RLQs. 

3. The constraints on Oqx evolution with redshift and luminosity can be greatly improved as the SDSS 
AGN sample increases in size. Ultimately it should be possible to perform analyses like those in VBS 
with ~ 3, 400 SDSS RQQs. Further X-ray studies of faint, moderate-luminosity AGN at z > 4 will 
also be very useful for breaking the degeneracy between redshift and luminosity. 

4. The current X-ray spectral constraints on z > 4 AGN are admittedly limited; these need to be 
improved. In addition to the underlying power law, low-redshift AGN show X-ray spectral complexity 
including neutral and ionized absorption, iron Ka lines, reflection continua, and soft excesses. The 
X-ray constraints on neutral absorption are already becoming interesting, but constraints on other 
types of spectral complexity are weak. With 40-100 ks exposures, XMM-Newton can obtain moderate- 
quality (~ 1000-2000 count) ~ 1-50 keV rest-frame spectra of luminous z > 4 quasars with 0.5-2 keV 
fluxes of >1.5 X 10^^^ erg cm~^ s~^ (see Figure lb). Efficient X-ray spectroscopy of the majority of 
the z > 4 AGN population, however, will require future missions such as Constellation-X , XEUS, and 
Generation-X (see Vignali et al. 2001, 2003a for relevant simulations). 

5. Systematic variability studies of z > 4 AGN would allow the tentative claim that quasar X-ray vari- 
ability increases with redshift (Manners et al. 2002) to be tested in extremis. 

6. Further X-ray imaging studies of z > 4 AGN will ultimately help to pin down the incidence of X-ray jets 
and gravitational lensing. Constraints will also be placed upon the amount of dust in the intergalactic 
medium and AGN clustering in large-scale structures. 
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